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Global attractivity and positive almost periodic 
solution of a discrete multispecies Gilpin-Ayala 
competition system with feedback control 
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Abstract — This paper discusses a discrete multispecies 
Gilpin-Ayala competition system. We first achieve the 
permanence of the system. Assume that the coefficients in the 
system are almost periodic sequences, we obtain the sufficient 
conditions for the existence of a unique almost periodic solution 
which is globally attractive. One example together with 
numerical simulation indicates the feasibility of the main results. 

Index Terms — Permanence , Almost periodic solution , 
Discrete, Gilpin-Ayala competition system, Global attractivity 

I. Introduction 

As we all known, investigating the almost periodic 
solutions of discrete population dynamics model with 
feedback control has more extensively practical application 
value(see [1-9] and the references cited therein). Wang [2] 
considered a nonlinear single species discrete with feedback 
control and obtain some sufficient conditions which assure 
the unique existence and global attractivity of almost positive 
periodic solution. Niu and Chen [5] studied a discrete 
Lotka-Volterra competitive system with feedback control and 
obtain the existence and uniqueness of the almost periodic 
solution which is uniformly asymptotically stable. 

In this paper, we investigate the dynamic behavior of the 
following discrete multispecies Gilpin-Ayala competition 
model with feedback controls 

Xi(k + 1) = Xi(k) exp ^ <Hj(k)(xj(k)) 0ii - ei (k)ui (k) j , 

< ( 1 - 1 ) 

n 

A Ui(k) - -fj (k)uj (k) + V 9ij ( k)xj (fc), i = 1,2, ,n, 

, j=i 

where {«,/£)}, {b,(k)}, {e,(k)}, {f(k ) ) and {gy(fc)} are 
bounded nonnegative almost periodic sequences such that 

0 < a l {j < a {j (k) < a£, 0 <H<b i (k)<b? 1 0 < e\ < e^k) < ej*, 

0<fl< Mk ) </“<!, 0 < < 9ij (k) < gf jt (1-2) 

hj= 1 . 2 , • • • , n, k^Z. 

Denote as Z and Z + the set of integers and the set of 
nonnegative integers, respectively. For any bounded sequence 
{g(/0} defined on Z, define 

g u = sup g(n),g l = info(ra). 
n£Z n € Z 

From the point of view of biology, in the sequel, we 
assume that x(0) = (xi(0), x 2 (0), • • • , x n (0), ^i(O), u 2 (0), • • • , 
u n ( 0)) > 0. Then it is easy to see that, for given x(0) > 0, the 
system (1.1) has a positive sequence solution x(k) = (xi ik), 
x 2 (k), • • • , x n (k))(k ^ Z + ) passing through x(0). 

With the stimulation from the works [10-13], the main 


purpose of this paper is to obtain a set of sufficient conditions 
to ensure the existence of a unique globally attractive positive 
almost periodic solution of system (1.1) with initial condition 
(1.3). 

The remaining part of this paper is organized as follows: In 
Section 2, we will introduce some definitions and several 
useful lemmas. In the next section, we establish the 
permanence of system (1.1). Then, in Section 4, we establish 
sufficient conditions to ensure the existence of a unique 
positive almost periodic solution which is globally attractive. 
The main results are illustrated by an example with numerical 
simulation in Section 5. Finally, the conclusion ends with 
brief remarks in the last section. 

II. Preliminaries 

First, we give the definitions of the terminologies involved. 

Definition 2.1([14]) A sequence x: Z^R is called an 
almost periodic sequence if the e-translation set of x 

E{s,x] = {t € Z :| x(n -\ - r) — x(n) |< £,Vn £ Z } 

is a relatively dense set in Z for all s > 0; that is, for any given 
£ > 0, there exists an integer /(e) > 0 such that each interval of 
length /(e) contains an integer x^E{e, x} with 

x(n -hr) — x(7?) |< s y Vn 6 Z. 

x is called an e-translation number of x(n). 

Denition 2.2( [15]) A sequence x : Z + — > R is called an 
asymptotically almost periodic sequence if 

x(n) = pin) + q(n), \/n £ Z + , 

where pin) is an almost periodic sequence and qin)^> 0, 

72— >00. 

Denition 2.3( [16]) A solution (xi(^), x 2 ik), • • • , x n ik)) of 
system (1.1) is said to be globally attractive if for any other 
solution (x i (&), x 2 ik), • • • , x n (&)) of system (1.1), we have 

Inn ( x * ik) — Xi(k)) = 0, i = 1, 2, • • • , n. 

k— M-oo 

Now, we present some results which will play an 
important role in the proof of the main result. 

Lemma 2.1([17]) If fx(zz) } is an almost periodic sequence, 
then fx(zz) } is bounded. 

Lemma 2.2([18]) fx(zz) } is an almost periodic sequence if 
and only if, for any sequence m/CZ, there exists a 
subsequence {m ik } such that the sequence {x(t7 +m ik ) } 
converges uniformly for all n ^Z as >o o. Furthermore, the 
limit sequence is also an almost periodic sequence. 

Lemma 2.3( [15]) fx(/r) } is an asymptotically almost 
periodic sequence if and only if, for any sequence mjC Z 
satisfying m^O and n^— >oo as /— >oo there exists a subsequence 
{m ik } c {mi} such that the sequence {x(n + m ik )} converges 
uniformly for all n £ Z + as k— >oo. 
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Lemma 2.4([ 17]) Suppose that {pi(n)| and {p 2 (n) } are 
almost periodic real sequences. Then {px(n)+p 2 (n)} and 
{pi(n)p 2 (n)} are almost periodic; l/pi(n) is also almost 
periodic provided that pi(n) 7 D for all n ^Z. 

Moreover, if e > 0 is an arbitrary real number, then there 
exists a relatively dense set that is ^-almost periodic common 
to { p 1 (n) } and {p 2 (n)J. 

Lemma 2.5( [19]) Assume that sequence {x(n)} satisfies 
x(n) > 0 and 

x(n + 1) < x(n) exp{a(n) — b(n)x a (n)} (2.1) 

for n^N, where a( n) and b( n) are non-negative sequences 
bounded above and below by positive constants. Then 

limsupx(n) < ( — rr ] Q exp{a u — — (2.2) 
n— >+oc \Gtb J OL 

Lemma 2.6( [20]) Assume that sequence {x(n)} satisfies 
x(n -h 1) > x(n) exp{a(n) — b(n)x ct (n)} , n > No , 


limsup < x* 

71 — t ■)“ OO 


and x(N 0 )> 0, where a( n) and b( n) are non-negative sequences 
bounded above and below by positive constants and N 0 ^ N. 
Then 


lint inf x(n) > min 

n — »+oo 



1 

C* f l 
exp {a 


b u (x *) a }, 



(2.3) 


Lemma 2.7( [21]) Assume that A > 0 and y(0) > 1, and 
further suppose that 

y(n + 1 ) < Ay(n) + B(n), n = 1 , 2 , 3 , • • • . 

Then for any integer k < n, 

k - 1 

y(n) < A k y(n — k) + A l B(n — i — 1). 

i=0 

Especially, if A < 1 and B is bounded above with respect to M, 
then 


limsup t/ (7?) < 


M 

l — A 


Lemma 2.8( [21]) Assume that A > 0 and y(0) > 1, and 
further suppose that 

y(n + 1) > Ay(n) + B(n), n = 1,2,3, ••• . 

Then for any integer k < n, 

fc-i 

y(n) > A k y(n - k ) + ^ A 1 B(n — i— 1). 

i=Q 

Especially, if A < 1 and B is bounded below with respect to m, 
then 

lim inf t/(n) > — . 

71— >OC 1 — 


nii < lini inf Xi(k) < lim sup Xi(k) < A f*, 

k—>~k-oc A:— >-+oo 

Hi < lim inf Ui(k) < lim sup Ui(k) < i = 1,2, • • • , n. 

k ~ ^+°0 

Proof. Let (x^k), x 2 (k), • • • , x n (k), ui(k), u 2 (k), • • • , u n (k)) 
be any positive solution of system (1.1). From the first 
equation of system ( 1 . 1 ), it follows that 

Xi(k+ 1) < Xi(k) exp |6j(A) - aa(k)(xi(k)) 0ii 1. (3.2) 


Thus, as a direct corollary of Lemma 2.5, according to (3.2), 
one has 

1 

lim sup Xi(k) < (— ^-p) exp ~ X = > 0. (3.3) 

fc-H- oo Uii a ii l J 

For any small positive constant 8 > 0, from (3.3) it follows 
that there exists a positive constants Kx > 0 such that for all k > 
Kx and i = 1 , 2 , • • • , n, 

Xi(k) < Ali + £. (3.4) 


For k > K b from (3.4) and system (1.1), we have 

n 

Ui(k + 1) < (1 — fl)ui(k) + ) " gjj(k)(Mj + £"). (3.5) 

j = i 

Then, as a direct corollary of Lemma 2.7, according to (3.5), 
one has 


1 

lim sup Ui (k) < -y ^ g™ (M j As). 

fc-^+oo Ji . = 1 

Letting 8^0, it follows that 

i n 

limsup Ui(k) < —j gjjMj = N ( . 

fc— >+o o J{ 


(3.G) 


Thus, there exists a positive integer K 2 > K 1? we have for k > 
K 2 


Ui(k) < Nt + s. (^* 7 ) 

For k > K 2 , from (3.7) and system (1.1), we have 

Xi(k+l) > ^(fcJexpjftiW - OiiikXxiik)) 9 " -ei(fc)(JVi+e)J. (3.8) 

Assuming that b\— c u \N[ > 0, for any 8 > 0, there exists a 
positive integer K 3 > K 2 such that bi(k) -ei(k)(Ni + 8 ) > 0 for 
k > K 3 . Thus, as a direct corollary of Lemma 2.6, according to 
(3.8), one has 


where 


lim inf Xj (k) > min{m 2ie , 7?ij 2£ }. 

k — >■ -|-oo 


m il£ = 


M — e “(Ni + e) \ g H 


a ti 


m i2(r = m il£ exp{6' -ef(Ni + e) - a“ i (M i ) 0 “ }. 


Letting e — » 0, it follows that 

lim inf Xi(k) > - min{m tl , m^} = rrii > 0, (3.9) 

k—y + oo 2 


where 


III. Permanence 

In this section, we establish the permanence result for 
system ( 1 . 1 ). 

Proposition 3.1 Assume that the conditions (1.2) and (1.3) 
hold, furthermore, 

6* — e“ATj > 0, (3.1) 

then system ( 1 . 1 ) is permanent, that is, there exist positive 
constants m i? M i? n^ and N A (i = 1, 2, • • • , n) which are 
independent of the solutions of system ( 1 . 1 ), such that for any 
positive solution (xx(k), x 2 (k), • • • , x n (k), Ui(k), u 2 (k), • • • , 
u n (k)) of system ( 1 . 1 ), one has: 


l 


m i2 = m il exp {b\ - e“ N t - a“(M t )° 11 }. 

From (3.9), for any 8 > 0, there exists a positive integer K 4 > 
K 3 such that 

Xi(k) >nii-£ (3.10) 

for k > K 4 . 

From (3.10) and system (1.1), we have 

n 

ui(k + 1) > (1 - fi)ui{k ) + y^gij(k)(mj - e). (3.11) 

j= i 

Then, as a direct corollary of Lemma 2.8, according to (3.11), 
one has 
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lim inf u l (k) > 

k — ^-|-oo 



Letting 8—4), it follows that 
l n 

lim inf u { (k) > — ^ g^rrij = m > 0. 

rt ) I <30 / ' 

j = 1 


(3.12) 


Then, (3.3), (3.6), (3.9) and (3.12) show that system (1.1) is 
permanent. The proof is completed. 

Theorem 3.2 Assume that (1.2), (1.3) and (3.1) hold, then 
system (1.1) is permanent. 


We denote by Q the set of all solutions (x^k), x 2 (k), • • • , 
x n (k), Ui(k), u 2 (k), • • • , u n (k)) of system (1.1) satisfying mi < 
Xi(k) < Mi, ni < Ui(k) < Ni(i = 1, 2, • • • , n) for all k ^ Z + . 

Proposition 3.3 Assume that (1.2), (1.3) and (3.1) hold. 
Then 

Proof. By the almost periodicity of { aij(k) } , {bi(k) } , { ei(k) } , 
{ fi(k) } and { gij (k) } , there exists an integer valued sequence 
{5 p } with 5 p — ► +oo as p — ► +oo such that 

aij(k + Sp ) — y ciij(k), bi(k 4- $ p ) — Y bi(k) y 

Qij(k + S p ) — > gij(k), ei(k + S p ) -A ei(k), fi(k + S p ) — > fi{k), 

as p — »• +cxj. 

Let 8 be an arbitrary small positive number. It follows from 
Theorem 3.3 that there exists a positive integer N 0 such that 

771 j — £ < Xi(k) < Mi + £, Jli — £ < Ui(k) < Ni + £, k > No. 

Write x ip (k) = Xi(k + 8 P ) and u ip (k) = Ui(k + 5 P ) for k > N 0 - 
5 p and p = 1 , 2, • • • . For any positive integer q, it is easy to see 
that there exists a sequence {x ip (k) : p > q} such that the 
sequence x p (k) has a subsequence, denoted by {x ip (k)} again, 
converging on any finite interval of Z + as p — > oo. Thus we 
have a sequence { yi(k) } such that 

x ip (k) -> Ui(k), Uj p (k) -4 Vi(k) for k e Z + as p -> oo. 
This, combined with 

.r,(k + 1 4- <5 p ) = x,(k + <5 p )exp < b,(k 4- <5 p ) — ^ <i, j(k + 6 p )(xj(k 4- S p )) 0,J — e,(k 4- S p )ii,(k 4- S p ) >, 

' j=l J 

i 

n 

&Ui(k + 6 p ) = —fj{k 4- S p )uj(k + 6 P ) + V g^(k + S p )xj(k + 6 P ), i = 1,2. ••• ,n, 
k j=i 

gives us 

Vi(k+ 1) = yi(k)expibi(k) - ^ a ij (k)(y j (k)) 9ij - e i (fc)v i (/e)|, 

^ j = l > 

i 

n 

A Vi(k) = -fi{k)vi(k) + ^ gij(k)yj(k), i = 1,2, • • • ,n. 

k. i= l 

We can easily see that (yi(k), y 2 (k), • • • , y n (k), Vi(k), 
v 2 (k), • • • , v n (k)) is a solution of system (1.1) and n^-s < 
yi(k) < Mi + 8, ^ - 8 < Vi(k) < Ni + 8 for k £ Z + . Since 8 is an 
arbitrary small positive number, it follows that mi < y^k) < M i? 
^ < Vi(k) < Ni and hence we complete the proof. 

IV. Global attractiyity and Almost Periodic 

Solution 

The main results of this paper concern the global 
attractivity of almost periodic solution of system (1.1) with 
conditions (1.2), (1.3) and (3.1). 


Theorem 4.1 Assume that (1.2), (1.3), (3.1) and 


n n 

(HI) pi = max{|l - ^ l» 1 1 ~^ 9 ij a ij M j ij \} + e i < b 

3= 1 3= 1 

n 

<Pi = l - fj + gijMj <1, * = ,n, 

3= 1 

hold. Then any positive solution (Xi(k), x 2 (k), • • • , x n (k), 
Ui(k), u 2 (k), • • • , u n (k)) of system (1.1) is globally attractive. 

Proof. Assume that (pi(k), p 2 (k), • • • , p n (k), Vi(k), 
v 2 (k), • • • , v n (k)) is a solution of system (1.1) satisfying 
(1.2) and (1.3). Let 

Xi(k) = Pi(k) exp{qi(k)}, u { (k) = v^k) + Wi(k), i = 1,2,--- ,n. 
Since 

qi(k+l) = \nxi(k + 1) — lnpi(k + 1) 

n 

- \nxi(k) + bi(k) - a i j{k)(x j (k)) 6ij - e^tyu^k) 

j = l 

n 

- In Pi(k) - bj(k) + ^ a ij(k)(Pj(k)) Gij + e { (k)Vi(k) 

j=i 

n 

= Qi(k) ~ L a n( k )l( x i( k )) 0ii ~ (Pj( k )) e * j ] ~ e,( k )w t (k) 
j = 1 

n 

= Qi( k ) ~ Y a n( k )(Pj( k )) eij [(exp{(jj(fc)}) s,J - l] - ei(k)wi{k) 
j= 1 

n 

= Qi(k)( 1 - J2 [Pj(k) exp{Aj (k) qj (k) }}°'> ) - ei (k)wi(k), 

3 = 1 

where 2i(k) ^ (0, 1). 

Similarly, we get 

Wi(k + 1) = Ui(k + 1) — Vi(k -\- 1) 

n n 

= (1 - fi(k))ui(k) + Y,9, i (k)x i (k) - (1 - fi(k))v,(k) -Y.9ii(k)Pj(k) 

j= 1 J=1 

71 

= (1 - fi(k))wi(k ) + ^ 9ij(k)Pj(k)(exp{qj(k)} - 1) 

j= 1 

71 

= (1 - fi(k))wi(k) + '^gij(k)p j (k)q j {k)exp{£, j (k)q j (k)}, i = 1, 2, — ,n, 

j=i 

where £j(k) £ (0, 1). 

To complete the proof, it suffices to show that 

lim qi(k) = 0, lim Wi(k) = 0, i = l,2,---,n. (4-1) 

k — >-\-oc k — ^ -foo 

In view of (HI), we can choose 8 > 0 such that 

71 71 

Pi = max{;l - ~ I, I - V (-< ;* < 1- 

j= 1 j = 1 


<Pi — 1 — fi + + s) < 1, i 1,2,*** , ?! . 

3= 1 

Let 

p = max{pf,</>f}, 

then p < 1. According to Theorem 3.2, there exists a positive 
integer k 0 ^ Z + such that 

nrii — £ < Xi(k) < + 8 , nii — £ < Pi(k) < Mi + £, 

i = 1, 2, • • * ,n 
for k > k 0 . 

Notice that ^(k) £ [0, 1] implies that pi(k) exp {Xi(k)ui(k) } 
lies between pi(k) and Xi(k), Lj(k) £ [0, 1] implies that pj(k) 
exp{Lj(k)Uj(k)} lies between pj(k) and Xj(k). From (4.1), we 
get 

n 

I Qi{k + 1)| < max{|l - ^ Oija^mj - £)° u \, 

j = i 

n 

|1 - Y*«“«(^i+e)*"l}l«(fc)|4e?Mfc)|, 

j= 1 
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M*+ 1)1 < (1- /')M*)I + y>“ j (M i +e)|« i (fc)|, 

j = 1 


i = 1, 2, • • • , n, for k > k 0 . 

In view of (4.3), we get 

max { max q t (k + 1 ) | , max | ( k + 1 ) | } < 

1 <i<n 


pmax{ max max 

KKn Ki<n 


W 


(*)!}, 


This implies 

max {max \qt(k)\ max K(fc)|} < 

l<i<n l<i<n 


k > k 0 . 


p k k ° max { max |^(A’o)|, max 1 1 l\ ( Ar 0 ) I } - 

l < Ci<n Ki<n 


Then (4.2) holds and we can obtain 
lim \x i{k) — pi(k) = 0, 

fc— H |-oo 


k > k 0 . 


lim \ui(k) — Vi(k) =0, i = 1,2, (4.2) 

k—y+oo 

Therefore, positive solution (x^k), x 2 (k), • • • , x n (k), Ui(k), 
u 2 (k), • • • , u n (k)) of system (1.1) is globally attractive. 


Theorem 4.2 Assume that (1.2), (1.3), (3.1) and (HI) hold. 
Then system (1.1) admits a unique almost periodic solution 
which is globally attractive. 

Proof. It follows from Proposition 3.1 that there exists a 
solution (xKk), x 2 (k), • • • , x n (k), ui(k), u 2 (k), • • • , u n (k)) of 
system (1.1) satisfying mi < X{(k) < M i? n{ < U{(k) < N i? k £ Z + . 

Suppose that (x^k), x 2 (k), • • • , x n (k), Ui(k), u 2 (k), • • • , 
u n (k)) is any solution of system (1.1), then there exists an 
integer valued sequence {k p }, k p — > +co as p — > +oo, such that 
(x ; (k + k p ), x 2 (k + k p ), • • • , x n (k +k p ), Ui(k + k p ), u 2 (k + 

? f 

k p ), • • • , u n (k + k p )) is a solution of the following system 

x i(k + 1) = Xi(k) exp | bi(k + k p ) - Y c iij(k + k p )(xj(k)) 9ij - e*(fc + fc p )ui(fc) j, 

^ i = i ' 

< 

n 

Aui(k) = -fi(k + k p )iii(k) + ^ gij(k + k p )xj(k), i = 1,2,**- ,n, 

, j =1 

From above discussion and Theorem 3.2, we have that not 

only (x^k+kp), x 2 (k+k p ), • • • , x n (k+k p ), Ui(k+k p ), u 2 (k + 

! » ? 

k p ), • • • , u n (k + k p )) but also (Ax^k + k p ),Ax 2 (k + 
k p ), • • • ,Ax n (k + k p ),Aui(k + k p ),Au 2 (k +k p ), • • • ,Au n (k + k p )) 

» I 

are uniformly bounded, thus (x^k + k p ), x 2 (k + k p ), • • • , x n (k 

? ! ! ? 

+ k p ), Ui(k + k p ), u 2 (k +k p ), • • • , u n (k + k p )) are uniformly 
bounded and equi-continuous. By Ascoli’s theorem[22], there 
exists a uniformly convergent subsequence (x^k+kp), 
x 2 (k+k p ), • • • , x n (k+k p ), Ui(k+k p ), u 2 (k+k p ), • • • , u n (k+k p )) 

£(xi(k + k p ), x 2 (k + k p ), • • • , x n (k + k p ), u^k + k p ), u 2 (k + 

» ? 

k p ), • • • , u n (k + k p )) such that for any 8 > 0, there exists a 
ko(s) > 0 with the property that if m, n > k 0 (s) then 

| X{(k + k m ) ~ Xi(k -b k n )\ < e, | Ui(k + k 7n ) - Ui(k + k n )\ < £, 

which shows from Lemma 2.2 that (Xi(k+k n ), x 2 (k+k n ), • • • , 
x n (k+k n ), Ui(k+k n ), u 2 (k+k n ), • • • , u n (k+k n )) is asymptotically 
almost periodic sequence, then (x^k + k n ), x 2 (k + k n ), • • • , 
x n (k + k n ), Ui(k + k n ), u 2 (k +k n ), • • • , u n (k+k n )) are the sum of 
an almost periodic sequence (pi(k+k n ), p 2 (k+k n ), • • • , 
p n (k+k n ), vi(k+k n ), v 2 (k + k n ), • • • , v n (k + k n )) and a sequence 
(qi(k + k n ), q 2 (k + k n ), • • • , q n (k + k n ),Wi(k + k n ),w 2 (k 
+k n ), • • • ,w n (k + k n )) defined on Z, such that 

Xi(k + k n ) = pi(k + k n ) + qi(k + k n ), 

Ui(k + k n ) = Vi(k + k n ) + Wi(k + k n ), k £ Z, 


where 

lim pi(k + k n ) = Pi{k ), lim Vi(k -b k n ) = Vi(k ), 

n— H-oo n— >+oo 

lim qj(k + k n ) =0, lim Wi(k + k n ) = 0, 

n—y+oo n — »-+oo 

(pi(k)} and { v A (k) } are almost periodic sequences, i = 1, 
2, • • • , n. It means that 

lim x i(k -b k n ) = Pi(k ), lim Ui(k -b k n ) = Vi(k). 

n. — y “l - oo 7i — y “I - oo 

In the following we show that { (pi(k), p 2 (k), • • • , p n (k), 
Vi(k), v 2 (k), • • • , v n (k))} is an almost periodic solution of 
system (1.1). 

From the properties of an almost periodic sequence, there 
exists an integer valued sequence {5 P }, 5 P — ► +oo as p — > +oo, 
such that 

bi(k ~b Sp ) ^ bi(k), ciij^k ~b Sp) ^ kiij^k), 

e i{k + S p ) —y ei(k ), fi(k + 5 p ) —y /;(fc), 

9ij(k + S p ) -a gij(k), as p -y -boo. 

It is easy to know that 

Xi(k -b S p ) -a Pi(k), Ui(k -b S p ) -a Vi (k) 
then we have 

Pi(k + 1) = lim Xi(k + 1 + (L) 

p—y oo 

= lim Xi(k + 8 P ) exp i bt(k + S p ) - aij(k + S p )(xj(k + S P )) 0U - e.i(k + S p )ui(k + S p ) i 

p—y oo I L ' 

v j — 1 J 

= pi(k) exp | bi(k) - a ij( k )(Pi( k )) e ' 3 ~ 

l j= i J 

Vi(k + 1) = lim Ui(k + 1 + S p ) 

p—y oo 

— p lim^ / [1 — fi(k + S p )]ui(k + 5 P ) + ^ ^ g t j (k + S p )xj(k + S p ) \ 

{ j'= i ' 

n 

= [1 ~ fi(k)]vi{k) + 'Yg ij (k)p j {k). 

3 = 1 

This prove that p(k) = {(pi(k), p 2 (k), • • • , p n (k), Vi(k), 
v 2 (k), • • • , v n (k))} satisfied system (1.1), and p(k) is a 
positive almost periodic solution of system (1.1). 

Now, we show that there is only one positive almost 
periodic solution of system (1.1). For any two positive almost 
periodic solutions (pi(k), p 2 (k), • • • , p n (k), Vi(k), v 2 (k), • • • , 
v n (k)) and (z^k), z 2 (k), • • • , z n (k), /{(k), / 2 (k), • • • , / n (k)) of 
system (1.1), we claim that pi(k) = Z{(k), V{(k) = /{( k) (i = 1, 
2, • • • , n) for all k £ Z + . 

jj* 

Otherwise there must be at least one positive integer K £ 
Z + such that pi(K ) ± Z{(K ) or Vj(K ) ± /j(K*) for a certain 
positive integer i or j, i.e., = lpi(K )~Z{(K )| > 0 or Q 2 = 

lvj(K ) _ /j(K )l > 0. So we can easily know that 

IF = lim pi(K*+6 p )- lim Zi(K*+6 p )\ = 

p—y+oo p-y- l-oo 

lim | pi (I\* -b S p )~ Zi (K* -b S p ) \ = lim \pi(k) - Zi(k) \ > 0, 

p — >-+oo k-y+oo 

or 


II 2 = 

lim 

p-y - f oo 


lim v j ( K * -b Sp ) — lim 

p—y+oo p— >•+ oo 

\ v j (K* + S p ) — lj (K* + i5 p ) | = 


lj{K* +d'p)| = 

= lim \vj(k)-lj(k)\ >0, 

k—y+oo 


which is a contradiction to (4.4). Thus pi(k) = Z{(k), V{(k) = 
/i(k)(i = 1, 2, • • • , n) hold for Vk^Z + . Therefore, system (1.1) 
admits a unique almost periodic solution which is globally 
attractive. This completes the proof of Theorem 4.2. 


V. Example and NUMERICAL SIMULATION 

In this section, we give the following example to check the 
feasibility of our result. 
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Example Consider the following almost periodic discrete 
Gilpin- Ayala competition system with feedback controls 

x\(k + 1) = xj(fc)exp 1 1.1 — 0.022 sin(>/3A) - (1.05 + 0.013sin(>/5A))x,(Jfc) 

- (0.025 - 0.001 cos(v/2 *))x 2 (*) - (0.02 + 0.0015cos(v'3*))x 3 (Jfc) ~ (0.025 - 0.002cos(v / 3Ar))u,(A*)|, 

x 2 (* + 1) = ar 2 (Jt)exp 1 1.15 - 0.025sin(v^*r)- (1085 + 0.015sin(v^3fc))x a (O 

- (0.025 + 0.003 <•<*)( ))x,( A- ) - (0.025 — 0.002 c-on(v / 2A-))x 3 (A-) - (0.025 + 0.004 sin(\/2Jfc))ti 2 (A) j. 

x 3 (A + 1) = x 3 (A-)exp 1 1.25 - 0.03sin(\/5Jb) - (1.1 - 0.024si»(>/2A:))x 3 (A) (5.1) 

- (0.03 - 0.002cos(>/2Jfc))xi(*) - (0.028 + 0.0015 cos(>/3*))x 2 ( k) - (0.02 + 0.002 cos( v / 3A))u 3 (A)|, 

A m{k) = -(0.83 - 0.03 sin(v / 2fc))w 1 (k) + (0.015 + 0.005 sm(y/3k))xi{k) 

+(0.013 - 0.004 sin(V3k))x 2 (k) + (0.024 - 0.005 cos(v / 5fe))rc 3 (fc), 

A u 2 {k) = -(0.824 - 0.04sin(\/3fc))w2(fc) + (0.018 - 0.004si n(Vbk))xi(k) 
+(0.015 - 0.005 cos{V2k))x 2 (k) + (0.014 + 0.004 sm{V2k))x 3 (k), 


A u 3 {k) = -(0.836 - 0.035 cos(VEk))u 3 (k) + (0.017 - 0.006 cos(\/2fc)):ci(&) 
+(0.013 - 0.005 sin( y/Sk))x 2 (k) + (0.014 + 0.005 cos(V2k))x 3 (k). 

By simple computation, we derive 

pi « 0.193, p-2 ~ 0.341, p3 w 0.298, 
tpi « 0.012, p2 ~ 0.075, (f2 ~ 0.218. 

It is easy to see that the conditions of Theorem 4.2 are verified. 
Therefore, system (5.1) has a unique positive almost periodic 
solution which is globally attractive. Our numerical 
simulations support our results(see Figure 1). 



0 10 20 30 40 50 60 70 

time k 
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Global attractivity and positive almost periodic solution of a discrete multispecies Gilpin -Ayala competition system 

with feedback control 


FIGURE 1: Dynamic behavior of positive almost periodic 
solution (xi(k), x 2 (k), x 3 (k), Ui(k), u 2 (k), u 3 (k)) of system (5.1) 
with the three initial conditions(0. 99, 0.97, 1.2, 0.071, 

0.068, 0.044), (0.97, 1.01, 1.1, 0.062, 0.064, 0.055) and (1.03, 
1.05,1.05,0.069,0.053,0.062) for k £ [1, 70], respectively. 

VI. Concluding Remarks 

In this paper, assuming that the coefficients in system 
(1.1) are bounded non-negative almost periodic sequences, 
we obtain the sufficient conditions for the existence of a 
unique almost periodic solution which is globally attractive. 
By comparative analysis, we find that when the coefficients in 
system (1.1) are almost periodic, the existence of a unique 
almost periodic solution of system (1.1) is determined by the 
global attractivity of system (1.1), which implies that there is 
no additional condition to add. 

Furthermore, for the almost periodic discrete 
multispecies Gilpin-Ayala competition system with 
timedelays and feedback controls, we would like to mention 
here the question of how to study the almost periodicity of the 
system and whether the existence of a unique almost periodic 
solution is determined by the global attractivity of the system 
or not. It is, in fact, a very challenging problem, and we leave 
it for our future work. 
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